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I.  INTRODUCTION 


Design  of  electron-beam  systems  in  linear  microwave  devices 
involves  three  main  elements:  beam  formation,  beam  focusing,  and  beam 
collection.  The  quality  of  the  electron  beim  formed  by  the  gun,  the 
focusing  system,  and  the  collection  system  each  play  a  major  role  in 
overall  tube  performance.  The  ability  to  put  these  three  elements 
together  and  achieve  high  tube  efficiency  is  of  utmost  importance  in 
microwave  tube  design. 

This  ever  present  demand  for  higher  tube  efficiency  has  led  to 
numerous  developments  in  the  areas  of  beam  focusing  and  collection. 
Improvements  in  beam  transmission  have  been  realized  through  newly 
developed  theories  that  established  optimum  entrance  conditions  neces¬ 
sary  for  proper  beam  focusing.*'*  With  the  use  of  an  electron  optics 
code  and  electron  beam  analyzers,  precise  positioning  of  the  beam  mini¬ 
mum  in  the  magnetic  field  region  is  accomplished  with  very  high  accu¬ 
racy.”*  Developments  in  the  area  of  beam  collection  have  also  con¬ 
tributed  greatly  in  better  tube  efficiency.  Ongoing  work  in  beam  col¬ 
lection  design  has  led  to  a  number  of  very  efficient  collection  schemes 
that  drastically  reduce  the  power  loss  arising  from  heat  dissipated  in 
the  collector.  The  most  significant  of  these  developments  is  multistage 
depressed  collectors,  which  are  reported  to  achieve  efficiencies  in 
excess  of  80  percent. 

These  developments  demonstrate  the  effects  of  proper  beam  coupling 
on  tube  performance.  They  also  point  out  the  fact  that  the  majority  of 
work  in  dc  beam  analysis  has  been  concentrated  in  the  collector  and  the 
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cathode-anode  region.  This  brings  out  the  tunnel  region  as  another 
avenue  in  the  continuous  effort  to  Improve  tube  performance.  Articles 
written  on  beam  behavior  along  the  focusing  structure  stress  its  impor¬ 
tance  in  spent  beam  analysis.  Accurate  description  of  dc  beam  behavior 
along  the  focusing  structure  will  not  only  result  in  vast  improvements 
in  spent  beam  models  used  in  computer  simulation  of  multistage  depressed 
collectors,  but  will  also  drastically  increase  the  accuracy  of  beam 
transmission  predictions.  It  is  therefore  essential  to  include  beam 
behavior  in  the  tunnel  region  in  any  beam  analysis.  This  is  the  moti¬ 
vating  factor  behind  this  study  —  the  need  to  investigate  beam  behavior 
along  the  focusing  structure. 

In  the  study  of  beam  behavior  along  the  length  of  the  tube,  a  need 
arises  for  a  means  of  measuring  beam  quality.  In  the  accelerator  com¬ 
munity,  this  is  accomplished  through  the  concept  of  emittance  or  phase 
space  area.  Standard  emittance,  formally  defined  as  the  area  of  the 
beam  phase  space  ellipse  divided  by  ir,  has  been  used  by  accelerator 
physicists  in  their  study  of  emittance  growth  found  to  occur  in  particle 
accelerators . 

This  concept  has  been  carried  over  in  the  microwave  tube  industry 
through  a  similar  parameter  introduced  by  Dr.  Richard  True  in  his  design 
of  gridded  guns.®  He  came  up  with  a  quantity  called  "tunnel  emittance," 
which  is  essentially  half  the  standard  emittance  normalized  with  respect 
to  the  beam  tunnel  radius.  In  his  recent  paper  on  tunnel  emittance,*  he 
used  this  parameter  to  intercouple  problems  generally  encountered  in 
electron  beam  systems  of  linear  beam  devices.  Based  on  his  observations 
that  indicate  occurrence  of  emittance  growth  in  PPM  focused  TWT's,  he 
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suggested  that  this  growth  could  explain  two  phenomena  frequently 


encountered  in  TWT's; 

1.  Degradation  in  beam  transmission  with  tube  length. 

2.  The  increased  rate  of  beam  spreading  of  a  dc  beam  after  pass¬ 
ing  through  a  PPM  structure. 

His  paper  emphasized,  among  other  things,  the  importance  of  tunnel 
emittance  in  the  design  and  analysis  of  beam  optical  systems. 

In  order  that  occurrence  of  emittance  growth  could  be  properly 
verified,  it  is  critical  to  have  the  ability  to  trace  the  electron  beam 
from  the  cathode  to  the  collector.  In  confined  flow  focused  tubes,  this 
could  be  easily  accomplished  by  optic  codes  with  the  aid  of  a  very  large 
computer.  The  same  could  not  be  said  in  the  case  of  PPM  focused 
tubes.  A  stability  problem  arising  In  the  solution  process  makes  it 
necessary  to  break  up  the  beam  and  solve  the  problem  In  segments.  This 
segmentation  process  poses  several  problems  with  regards  to  the  manner 
the  beam  Is  segmented  and  the  boundary  conditions  that  must  be  enforced 
to  ensure  a  consistent  solution.  No  criteria  have  yet  been  established 
with  regard  to  this  segmentation  process. 

Today  there  are  several  optics  codes  available  that  could  be  used 
in  ray  tracing.  One  such  code  is  the  deformable  mesh  electron  optics 
code  of  Dr.  True  (DEMEOS).*  It  has  found  use  In  both  electron  gun  and 
multistage  depressed  collector  design.  Numerical  results  of  electron 
gun  simulations  using  this  code  show  very  close  agreement  with  experi¬ 
mental  data  obtained  from  the  Litton  electron-beam  analyzer  system.** 
Numerous  computer  simulations  have  been  made  with  this  code  to  trace  the 
beam  far  downstream,  but  no  known  effort  has  been  made  to  study  the 
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accuracy  of  the  results  and  to  analyze  the  ray  tracing  method.  It  is 


imperative  that  the  ray  tracing  method  we  use  be  very  accurate,  since 
tunnel  emittance  is  an  extremely  sensitive  parameter.  It  is  the  goal  of 
this  project  to  study  and  develop  an  accurate  method  of  downstream  ray 
tracing  that  will  allow  verification  of  tunnel  emittance  growth  and 
enable  dc  beam  analysis  along  the  entire  length  of  the  tube.  The  DEMEOS 
electron  optics  code  will  be  studied  and  modified  to  perform  the  beam 
segmentation  necessary  in  PPM-focused  tube  simulations. 


II.  COMPUTER  MODELING 


This  section  presents  the  modeling  process  undertaken  in  using  the 
deformable  mesh  electron  optics  code  (DEMEOS).  Brief  discussions  of 
each  step  are  included  to  give  the  reader  a  general  idea  of  how  the 
model  was  obtained.  The  reader  should  refer  to  reference  3  for  a 
detailed  description  of  the  process. 

A.  Electron  Gun  Selection 

The  need  to  experimentally  verify  results  of  this  study  places 
some  restrictions  on  the  choice  of  gun  to  be  used  In  the  computer  simu¬ 
lation.  The  gun  must  provide  a  small  and  very  laminar  beam  to  be  able 
to  achieve  the  highest  beam  transmission  possible.  The  choice  of  design 
Is  also  limited  to  Litton  production  guns  to  assure  parts  availability, 
minimize  design  time,  and  lower  production  cost.  These  considerations 
led  to  the  choice  of  the  L5637  gun  for  this  study. 

B.  Gun  Modeling 

The  first  step  in  gun  modeling  using  the  DEMEOS  code  is  to  trace 
out  the  logical  diagram.  This  involves  outlining  the  gun’s  boundary 
contours  (shown  in  Fig.  1 )  In  a  regular  mesh  denoted  "logical"  space. 
Each  of  the  corner  points  In  the  logical  space  corresponds  to  a  real 
geometrical  position  and  are  associated  with  a  given  potential.  Figure 
2  shows  the  logical  diagram  for  the  L5637  gun.  Transformation  from  this 
logical  space  to  real  space  is  automatically  accomplished  by  the  code 
resulting  In  an  Irregular-triangle  relaxation  mesh  shown  In  Fig.  3.  It 
Is  on  this  relaxation  mesh  that  the  finite-difference  approximation  of 
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Fig.  1.  Electrode  contours  for  the  L5637  gun. 
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Logical  diagram  for  the  L5637  gun. 


Poisson’s  equation  is  solved.  Figure  shows  the  solution  for  this 
particular  gun. 

Once  the  model  of  the  gun  is  set  up  in  the  computer,  several 
simulation  runs  at  various  cathode-to-anode  distances  with  different 
grid  voltages  are  made.  From  these  runs,  hairpin  curves  shown  in  Fig.  5 
are  obtained.  Through  these  curves,  the  optimum  operating  perveance, 
cathode-to-anode  distance,  and  grid  voltage  for  the  gun  design  are 
determined. ® 

C .  Electron  Beam  Focusing 

This  gun  design  must  now  be  matched  to  a  suitable  focusing  struc¬ 
ture.  The  appropriate  magnetic  field  needed  to  focus  the  beam  formed  by 
this  gun  must  be  determined.  The  peak  magnetic  field  value  is  obtained 
using  True’s  code  GUNPPMC,  a  PC  program  based  on  his  work.®  The  print¬ 
out  of  the  output,  shown  in  Fig.  6,  shows  that  ^  ,^25  kC  rms  will  be 
sufficient  to  properly  focus  this  beam.  This  leads  to  the  choice  of  the 
L5638  tube  body  for  this  particular  gun  design.  Figure  7  shows  the 
axial  magnetic  field  profile  of  the  L5638  showing  a  peak  value  of  2.8 
kC.  The  values  used  in  the  simulation  were  obtained  from  the  actual 
values  of  the  L5638  magnetic  field  scaled  to  get  the  desired  peak  B- 
field  values.  The  first  single-peak  field  was  adjusted  and  positioned 
according  to  the  method  described  in  reference  1 .  Figure  8  shows  the 
resulting  magnetic  field  profile. 

Note  that  the  nonthermal  model  was  used  in  all  of  the  simulations 
in  this  study.  Simulations  were  made  using  thermal  and  nonthermal 
models  to  determine  the  agreement  between  the  resulting  solutions.  The 
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R95/A 


CiUNPPMC  (GUNPPMX):  ProgrAn  After  TrAtis.  E.D.  1984|  1987 
Vers.  1.3,  VAX,  BasIc,  R.  True,  10-23-87 


Units;  enter  1  for  cm,  2.54  for  inches?  2.54 
CAse  NAme  *?  L-5637  GUN  FOR  RAY  TRACING  STUDY 
enter  disk  rc,  VO  (kv),  i*mu?  .2,26, .7548 
enter  r95/A,  a  (ckt  i.r.J,  period  L 

(note:  pulse  r95/A  *  .6  opt,  .65  ax;  cw  .45  opt,  .5  ax) 

?  . 5463, .076, .524 

iO  =  3.1644  rb»  .415188E-01  rc/r95  *  4.81709 

caIc  th  *  24.5979  deg 

enter  actuAl  th  or  0  to  bypAss?  24.2 

Actual  th  X  24.2  deg  deltA  th  x  ,397911  deg 

CAth  sph  rad  a  .487896  cath  area  x  ,847993  ca*2 

ungr:  jO  x  3.4037  ja  x  3.73164  je  x  4.09117  a/ca*2 
beam  j  x  90.5701  a/cm*2 

enter  scrn  factor  (diode=0,  sh.g.x.i7  (typ),  i.gx.i  (typ))?  .17 
enter  cath  usage  factor  (sh.g.x.7  (typ),  i.g.  and  u.l.t.xl  (typ), 
c . g. x , 85  ( typ ) )?  .7 

gr:  jO  X  5.85835  ja  x  6.42279  je  x  7.0416  a/ca‘2 

calc  signa  *  .0363  (note:  calc  sigaa  froa  actual  th  if  entered) 
enter  actual  sigaa  (e.g.  froa  analyzer)  or  0  to  bypass?  .033264 

enter  T(deg  C)  to  calc  thermal  sigaa  or  0  to  bypass?  1130 
thermal  sigma  x  .734543E-02 

enter  factor  F  (1.65  default)  (note:  for  axial  CCTWT  B  vals  use  Fxl.3)?  1.3 
B*/Bbr  *  1.29482 

Bbr  X  1.10059  B4 ( ras )  x  1.42507  kgauss 

B+  pk  s.p.  X  2.01535  B*  pk  d.p.  x  2.37512  (note:  B^  pk  d . p. xB* ( ras )/ . 6 ) 

focfl  X  2.37017  lamda_p/L  x  3.29157 

(note:  pulse  lamda_p/L>2 . 7  and  cw  laoda_p/L>3 . 2 ) 

Vs.b.(kv)  X  3.04469  Vs.b./VO  =  .117104 
T-emittance  =  .181721E-01 


Fig,  6,  Printout  of  output  of  GUNPPMC  using  parameters 
from  DEMEOS  thermal  insulation. 
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used 


nonthermal  models  yielded  beam  filling  factors  and  t-emittance  values 
that  were  very  close  to  those  obtained  using  thermal  models.  The  effect 
of  thermal  velocities  in  these  emittance  growth  studies  are  minimal 
compared  to  the  effects  of  the  grid  on  the  value  of  sigma.  This  close 
agreement  between  the  thermal  and  nonthermal  simulations  and  the 
increased  computer  time  Involved  in  thermal  simulations  led  to  the 
choice  of  the  nonthermal  model  as  the  vehicle  in  our  investigation. 
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•  III.  SIMULATION 


This  section  discusses  the  process  of  ray  tracing  using  the  DEMEOS 
code.  It  describes  the  current  method  of  simulating  dc  beam  behavior 
from  the  gun  to  the  collector.  A  simulation  of  a  PPM  focused  tube  is 
presented  to  demonstrate  the  problems  encountered  in  the  use  of  this 
method. 

A .  DEMEOS  Simulation 

As  mentioned  before,  simulation  of  the  entire  region  of  a  PPM 
focused  tube  could  not  be  performed  without  encountering  instability  in 
the  solution  process.  This  Instability  is  avoided  by  limiting  the 
length  of  the  region  Included  in  the  simulation.  This  led  to  the  pres¬ 
ent  method  of  dividing  the  entire  region,  from  the  cathode  to  the  col¬ 
lector,  into  segments  and  solving  them  separately. 

The  first  segment  in  the  simulation  includes  the  gun  and  a  portion 
of  the  tunnel.  Poisson’s  equation  is  solved  in  this  region,  and  ray 
data  at  the  end  of  the  tunnel  are  saved  and  used  in  the  next  region. 
The  segments  following  the  gun  region  are  continuations  of  the  tunnel. 
Continuation  runs  are  made  until  the  collector  region  is  reached.  The 
length  of  the  segments  in  these  continuation  runs  is  dictated  by  the 
following  restrictions: 

1  .  The  length  must  be  such  that  the  solution  will  be  conver¬ 
gent.  It  is  desirable  to  make  the  segment  as  long  as  possible 
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to  speed  up  the  solution  process.  On  the  other  hand,  increas¬ 
ing  the  length  of  the  segment  will  require  a  much  longer 
convergence  time. 

2.  Simulation  runs  must  end  and  begin  at  the  location  of  zero 
magnetic  field.  This  restriction  arises  because  DEMEOS  does 
not  take  into  account  any  rotational  motion  of  the  beam  at  the 
beginning  of  each  simulation  run.  This  was  a  deliberate 
effort  to  avoid  the  additional  complexity  that  the  beam's 
rotational  motion  will  bring  to  the  code.  Knowing  that  the 
beam's  rotational  motion  is  zero  at  the  location  of  zero 
magnetic  field,  the  beginning  and  ending  of  any  simulation  at 
the  location  of  zero  magnetic  field  Justifies  neglect  of  this 
rotational  motion.  This  limitation  does  not  in  any  way  affect 
the  accuracy  of  the  result. 

B.  Simulation  with  Beam  Segmentation 

A  simulation  run  using  the  present  segmentation  method  (with  more 
segments  than  is  customary)  was  made  and  the  resulting  beam  profile  is 
shown  in  Fig.  9.  The  length  of  the  first  segment  was  adjusted  to 
include  the  first  single-peak  magnetic  field.  Termination  was  made  to 
coincide  with  the  location  of  zero  magnetic  field.  The  succeeding 
segments  or  continuation  runs  were  all  made  to  each  include  one  full 
period  of  the  double-peak  magnetic  field. 

The  beam  profile  of  Fig.  9  shows  an  electron  beam  that  is  getting 
more  laminar  as  it  moves  away  from  the  cathode.  The  plots  of  the  tunnel 
emittance  and  beam  radius  along  the  length  of  the  tunnel,  shown  in  Figs. 
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Fig.  9.  Electron  trajectories  using  the  present  ray  tracing  method  with  beam  segmentation. 


10  and  11,  show  the  beam  decreasing  in  size  and  exhibiting  negative 
emittance  growth.  The  electrons  are  behaving  as  if  they  are  experienc¬ 
ing  a  cooling  down  process  as  they  flow  away  from  the  cathode.  To 
verify  the  validity  of  this  solution,  a  second  simulation  run  was  per¬ 
formed  in  a  different  fashion. 

C.  Simulation  without  Beam  Segmentation 

As  noted  before,  no  work  has  been  made  to  verify  the  accuracy  of 
the  present  ray  tracing  method.  On  the  other  hand,  results  of  simula¬ 
tion  runs  without  beam  segmentation  have  been  verified  to  be  accurate 
through  the  use  of  electron  beam  analyzers.  This  provided  a  means  by 
which  the  solution  of  the  first  simulation  run  could  be  verified.  By 
simulating  the  entire  length  of  the  first  simulation  without  breaking  up 
the  beam,  a  reliable  solution  is  obtained.  To  avoid  the  problem  of 
instability,  the  length  of  the  region  included  in  the  simulation  was 
adjusted  until  convergence  of  the  solution  was  reached. 

Figures  12,  13f  and  1^1  show  the  beam  profile  resulting  from  this 
simulation  and  the  corresponding  plots  of  tunnel  emittance  and  beam 
radius  along  the  tunnel.  The  figures  show  an  electron  beam  behaving 
completely  different  from  the  beam  of  Fig.  9.  The  beam  is  very  nonlami- 
nar  and  the  smoothing  observed  before  has  disappeared.  The  difference 
between  the  two  solutions  is  emphasized  by  superimposing  the  curves  as 
shown  in  Figs.  15  and  16.  Since  solutions  of  other  beam  optics  problems 
obtained  in  a  similar  process  have  been  proven  to  be  accurate,  the  beam 
profile  of  Fig.  12  is  the  true  picture  of  electron  behavior  inside  the 
tunnel . 
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Fig.  11.  Beam  filling  factor  versus  axial 
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emittance  versus  axial  distance  of  the  beam  profile  of  Fig.  12 


Fig.  Beam  filling  factor  versus  axial 


emittance  versus  axial  distance  for  the  two  electron  trajectories  of  Figs.  9  and  12 


Fig.  16.  Beam  filling  factor  versus  axial  distance  for  the 
two  electron  trajectories  of  Figs.  9  and  12. 


D.  Beam  Segmentation  Problem 


These  two  very  different  beam  profiles  for  the  same  electron 
optics  problem  clearly  manifest  the  existence  of  a  problem  In  the  pres¬ 
ent  method  of  ray  tracing.  This  demonstrates  that  the  process  of  split¬ 
ting  the  region  Into  small  segments  and  solving  them  Individually  Is  not 
a  trivial  process.  Guidelines  on  proper  beam  segmentation  need  to  be 
established  to  make  this  ray  tracing  method  accurate  and  reliable. 
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IV.  STUDY  OF  THE  RAY  TRACING  METHOD 


The  previous  section  pointed  out  the  presence  of  a  problem  in  the 
present  ray  tracing  method.  The  next  task  is  to  determine  the  cause  of 
the  problem.  This  section  will  deal  with  these  matters.  The  present 
ray  tracing  method  will  be  investigated  and  the  parameters  in  the  simu¬ 
lation  process  will  be  studied  to  determine  their  effects  on  the  solu¬ 
tion.  Possible  sources  of  error  will  be  discussed,  and  the  role  they 
play  in  altering  the  solution  will  be  pointed  out. 

A.  Boundary  Specifications 

Determination  of  the  approximate  area  where  the  problem  initially 
surfaced  is  accomplished  through  careful  study  of  the  beam  profiles  of 
Figs.  9  and  12.  Close  examination  reveals  the  profiles  almost  identical 
in  the  first  region,  but  very  much  different  in  the  following  region 
(see  Figs.  17,  18,  and  19).  This  suggests  that  the  problem  lies  in  the 
extraction  and  injection  of  ray  data  from  one  simulation  run  to  another. 
Focus  of  attention  could  therefore  be  directed  to  these  processes. 

One  factor  that  plays  a  major  role  in  the  solution  process,  espe¬ 
cially  in  the  ray  data  extraction  and  Injection,  is  the  boundary  speci¬ 
fications.  It  is  at  the  boundaries  that  rays  are  extracted  and 
injected.  It  is  also  through  the  conditions  specified  at  the  boundaries 
that  the  finite-difference  approximations  of  the  Laplace's  and  Poisson's 
equations  are  solved.  It  is  crucial  that  the  proper  boundaries  be 
correctly  specified  to  ensure  consistency  in  the  solution. 
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vs  Z  (centimeters) 


The  three  main  boundaries  used  in  the  DEMEOS  code  are  the  B,  C, 
and  N  boundaries.  These  boundaries  are  shown  on  the  logical  diagram  of 
Fig.  2.  The  C  boundary,  normally  associated  with  the  cathode  surface, 
is  not  of  interest  to  us  in  this  study,  since  it  is  not  involved  in  data 
extraction  or  injection.  The  B  and  N,  however,  are  the  major  boundaries 
involved  in  the  segmentation  process. 

B  boundaries  are  boundaries  that  apply  to  regions  associated  with 
fixed  potentials.  It  is  synonymous  with  the  Dirichlet  boundary  wherein 
the  potential  function  is  known  on  the  boundary  surface.  Examples  of 
this  boundary  are  the  focus  electrode  contour,  the  anode,  and  the  tunnel 
wall.  N  or  Neuman  boundaries,  on  the  other  hand,  apply  to  regions  where 
the  normal  derivative  of  the  potential  function  is  known.  They  are 
shown  on  the  logical  diagram  defining  regions  where  symmetry  is  assumed 
to  exist.  Since  the  normal  derivative  is  also  an  unknown  in  the  prob¬ 
lem,  DEMEOS  uses  the  special  case  of  the  Neuman  boundary  condition, 
i.e.,  the  normal  derivataive  of  the  potential  function  is  set  to  zero. 
This  in  effect  forces  the  potential  to  remain  constant  across  the  bound¬ 
ary.  This  plays  a  very  significant  role  in  the  solution  process. 

B.  Neuman  Boundary  Effects 

To  determine  the  extent  of  the  effect  of  this  artificial  condition 
enforced  across  the  boundary,  a  study  was  made  in  which  a  very  laminar 
beam  with  a  uniform  current  density  profile  was  injected  inside  a  tunnel 
whose  boundaries  are  defined  according  to  the  logical  diagram  of  Fig. 
20a.  Simulation  was  made  under  the  following  conditions: 
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(a)  Tunnel  for  Case  I. 
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(b)  Tunnel  for  Case  II. 


Fig.  20.  Logical  diagram  of  the  tunnel  used 
in  the  Neuman  boundary  study 


Vo  26  kV 

Perveance  0.755  microperveance 

Tunnel  length  1  cm 

Tunnel  radius  0.16M28  cm 

B-fleld  start  10  kG 

The  beam  trajectory  was  controlled  by  varying  the  magnetic  field 
linearly  along  the  tunnel.  The  magnetic  field  was  adjusted  to  obtain 
beams  with  negataive,  positive,  and  zero  trajectories.  Two  cases  were 
simulated: 

Case  I  The  entire  length  of  the  tunnel  was  included  in  the 
simulation. 

Case  II  Only  one  half  of  the  tunnel  region  of  Case  I  was 
Included.  The  boundaries  were  defined  according  to  the 
logical  diagram  of  Fig.  20b.  All  parameters  were  kept 
the  same  as  Case  I,  except  for  the  length  of  the  region 
Included  in  the  simulation.  This  change  was  made  to 
determine  the  effect  of  a  Neuman  boundary  in  the  middle 
section  (z  »  0.5  cm)  of  the  region  of  Case  I.  Any 
difference  between  the  solution  of  this  case  and  the 
solution  of  Case  I  could  be  attributed  solely  to  this 
Neuman  boundary . 

Figure  21  shows  the  beam  profiles  resulting  from  the  simulation  of 
the  entire  tunnel.  It  is  noteworthy  to  mention  that  each  of  these 
profiles  resulted  in  a  unique  potential  distribution.  Comparison  of  the 
solutions  of  the  two  cases  is  made  by  plotting  the  potentials  at  a 
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constant  nadius’ along  the  tunnel  region.  Figures  20,  21,  and  22  show 
the  plots  of  these  potentials  for  the  three  potential  gradients. 

The  solid  curve  of  Fig.  22  is  a  portion  of  the  potential  distribu¬ 
tion  at  a  constant  radius  along  the  tunnel  of  Case  I.  It  shows  the 
potential  remaining  constant  at  a  radius  of  r  -  0.1^66  cm  along  the 
tunnel.  The  dashed  curve,  shown  coincidental  with  the  solid  curve, 
represents  the  potential  distribution  at  the  same  radius  along  the 
tunnel  of  Case  II.  The  agreement  between  these  two  solutions  demon¬ 
strates  that  insertion  of  the  Neuman  boundary  at  z  »  0.5  cm  did  not 
change  the  solution.  This  is  as  expected,  since  the  condition  that 
Neuman  boundary  enforces  in  Case  II  happens  to  be  the  condition 
occurring  at  that  location  in  Case  I.  Case  II  Is  just  a  simulation  of 
the  first  half  of  Case  I.  We  could  therefore  conclude  that  the  solution 
is  not  affected  by  the  special  Neuman  condition  in  cases  where  the 
potential  distribution  is  constant  across  the  boundary.  Figure  23  shows 
the  case  where  the  potential  is  increasing  along  the  tunnel.  The  two 
curves  are  coincidental  only  at  regions  away  from  the  boundary.  They 
begin  to  deviate  from  each  other  as  they  approach  the  location  of  the 
Neuman  boundary  at  z  *  0.5  cm.  The  same  behavior  Is  exhibited  in  Fig. 
2^,  the  case  where  the  potential  is  decreasing  along  the  tunnel.  The 
difference  between  the  two  solutions  can  be  explained  by  noting  that  the 
condition  enforced  by  the  Neuman  boundary  in  Case  II  (i.e.,  the  poten¬ 
tial  is  constant  across  z  -  0.5  cm)  is  different  from  the  condition 
occurring  at  that  same  location  in  Case  I.  The  potential  is  varying  at 
that  location  of  the  tunnel  in  Case  I.  The  Neuman  boundary  is  clearly 
having  an  effect  on  the  solution.  The  imposition  of  a  Neuman  boundary 


-  36  - 


Fig.  22.  Potential  distribution  at  a  fixed  radius  inside  the  tunnel 
of  Fig.  20  (constant  potential  across  the  boundary). 


Fig.  23.  Potential  distribution  at  a  fixed  radius  inside  the  tunnel 
of  Fig.  20  (increasing  potential  across  the  boundary). 


Fig.  24.  Potential  distribution  at  a  fixed  radius  inside  the  tunnel 
of  Fig.  20  (decreasing  potential  across  the  boundary). 


at  locations  where  the  potential  is  not  constant  will  therefore  result 
in  an  erroneous  solution.  This  is  a  very  significant  finding.  This 
could  explain  the  beam  alteration  observed  in  the  segmentation 
process.  Ray  data  injected  in  the  continuation  runs  are  being  extracted 
at  the  Neuman  boundary,  and  in  cases  where  the  potential  is  not  constant 
across  the  boundary  (which  often  is  the  case),  these  ray  data  are  com¬ 
pletely  erroneous. 

This  new  discovery  regarding  the  Neuman  boundary  places  us  in  a 
very  difficult  predicament.  It  refrains  us  from  using  the  special 
Neuman  boundary  condition,  which  will  make  it  virtually  impossible  to 
solve  any  beam  optics  problem.  A  compromise  has  to  be  found  that  will 
permit  use  of  the  special  Neuman  boundary  and  still  lead  to  a  valid 
solution.  The  plots  of  Figs.  23  and  24  give  an  indication  that  such  a 
compromise  is  realizable.  The  curves,  shown  coincidental  at  some  dis¬ 
tance  away  from  the  boundary,  suggest  that  a  portion  of  the  solution  is 
unaffected  by  the  error  occurring  at  the  boundary.  Erroneous  data  could 
therefore  be  Isolated  by  discarding  the  portion  of  the  solution  close  to 
the  boundary.  We  only  have  to  determine  how  much  of  the  solution  next 
to  the  boundary  is  erroneous.  This  could  be  accomplished  by  determining 
the  distance  where  the  two  solutions  begin  to  deviate  from  each  other. 

To  make  a  determination  as  to  how  far  from  the  boundary  the  error 
propagates,  a  beam  with  a  positive  trajectory  was  again  simulated  with 
the  tunnel  length  of  2  cm.  Simulation  runs  similar  to  Case  II  were  then 
made  with  lengths  of  0.5,  1.0,  and  1.5  cm.  The  resulting  potential 
distributions  were  then  compared  to  determine  the  distance  the  solutions 
of  the  shorter  runs  begin  to  vary  from  the  solution  of  the  2  cm  tunnel. 
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The  resulting  distances  are  presented  in  the  plot  of  Fig.  25.  The 
length  of  the  region  unaffected  by  the  Neuman  boundary  (denoted  by  L') 
is  shown  to  increase  as  the  simulation  segment  length  L  is  increased. 
At  a  simulation  length  of  1.5  cm,  the  error  propagated  only  0.225  cm 
from  the  Neuman  boundary.  This  means  that  85  percent  of  the  solution  is 
unaffected  by  the  Neuman  boundary.  This  very  important  result  is  the 
answer  to  our  present  dilemma.  It  suggests  that  the  error  propagating 
from  the  Neuman  boundary  is  directly  dependent  on  the  length  of  the 
simulation.  Increasing  the  length  of  simulation  runs  will  result  in  an 
increase  in  the  valid  or  usable  portion  of  the  solution  and  allows  us  to 
continue  using  the  special  Neuman  boundary  condition.  (We  must  be 
reminded  that  the  simulation  lengths  are  still  limited  by  the  insta¬ 
bility  problem  discussed  earlier.) 

To  summarize  the  results  obtained  from  this  Neuman  boundary  study, 
the  significant  findings  are  enumerated  as  follows: 

1.  The  use  of  the  special  Neuman  boundary  in  beam  segmentation  is 
one  cause  of  error  in  the  solution  process  and  should  there¬ 
fore  be  avoided  as  much  as  possible. 

2.  The  solution  adjacent  to  the  Neuman  boundary  is  erroneous  and 
should  be  discarded. 

3.  Error  propagation  could  be  minimized  by  increasing  the  length 


of  the  simulation 


1.00  1.50 

(CM) 


2.00 


as  a  function  of  simulation  length 


C .  Magnetic  Field  Effeota 


Another  factor  that  was  found  to  affect  the  result  of  the  simula¬ 
tion  was  the  magnetic  field  at  the  start  of  the  simulation.  Section  III 
emphasized  the  Importance  of  starting  all  simulations  at  the  zero 
magnetic  field  point.  Strict  adherence  to  this  restriction  is  very 
critical  in  the  solution  process.  When  simulation  runs  are  ended  at 
locations  other  than  where  B  Is  equal  to  zero,  the  beam  will  have  a 
rotation  that  will  not  be  accounted  for  at  the  start  of  the  following 
continuation  runs.  This  was  experienced  in  this  study. 

In  all  the  simulation  runs  made  in  this  study,  actual  magnetic 
field  values  were  used.  The  magnetic  field  analyzer  used  in  the 
measurement  was  set  to  print  out  the  field  values  every  5  mils.  This 
resulted  In  the  location  of  the  zero  magnetic  field  not  being  discretely 
defined.  Since  the  field  values  in  the  transition  from  +B  to  -B  were 
very  small  (from  +0.10176  G  at  2.80614  cm  to  -0.26288  G  at  2.8191  cm  for 
the  first  zero  crossing),  the  simulation  was  ended  at  the  location  where 
the  field  made  its  first  sign  change  (at  z  *■  2.8191  cm  for  the  first 
segment).  Assumption  was  made  that  the  B-field  value  of  -0.26288  G  was 
small  enough  that  any  error  occurring  from  it  will  be  negligible.  This 
assumption  had  a  very  drastic  effect  in  the  electron  trajectories.  This 
is  shown  in  Figs.  19a  and  19b.  Electrons  that  were  reaching  the  axis  in 
Fig.  19b  were  completely  diverted  from  their  paths  and  caused  to  flow 
towards  the  tunnel  wall.  This  resulted  in  a  decrease  in  beam  current  in 
each  continuation  run.  This  contributed  to  the  smoothing  in  beam  tra¬ 
jectory  and  decrease  in  beam  size  demonstrated  in  Figs.  9  and  10. 
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D.  R95  Values 


The  R95  values  printed  in  the  DEMEOS  output  were  also  found  to  be 
a  potential  source  of  error  without  proper  interpretation.  In  the  plots 
of  the  beam  filling  factor  versus  the  axial  distance,  R95  values  were 
obtained  directly  from  the  DEMEOS  output.  They  are  calculated  in  the 
code  based  on  the  beam  current  present  at  the  start  of  the  simulation. 
In  the  first  simulation  run,  this  beam  current  is  the  total  current 
originating  from  the  cathode  minus  the  amount  intercepted  by  the  grid. 
In  the  continuation  runs,  the  starting  current  is  the  total  beam  current 
from  the  cathode  minus  the  current  intercepted  by  the  grid  and  the 
tunnel  wall.  Since  R95  is  defined  as  the  radius  that  encloses  95  per¬ 
cent  of  the  total  beam  current  that  passes  through  the  grid,  the  R95 
values  of  continuation  runs  are  therefore  erroneous.  These  erroneous 
R95  values  in  the  continuation  runs  were  used  in  the  plots  of  Figs.  10 
and  11.  This  error  in  the  computation  of  the  R95  values  coupled  with 
the  increased  beam  interception  resulting  from  the  improper  initial  B- 
field  values  is  another  explanation  for  the  negative  emittance  growth 
and  decrease  in  beam  size  observed  in  the  beam  profile  of  Fig.  9. 


V.  MODIFICATIONS  ON  THE  BEAM  SEGMENTATION .METHOD 


This  section  deals  with  the  modification  necessary  to  make  the 
solution  of  the  present  ray  tracing  method  accurate  and  reliable.  It 
introduces  a  new  technique  of  data  extraction  and  injection  that  greatly 
reduces  the  effect  of  the  error  emanating  from  the  Neuman  boundary.  It 
also  briefly  mentions  the  corrections  needed  to  account  for  errors  due 
to  the  magnetic  field  and  the  R95  values  in  the  continuation  runs. 

A.  Modifications  and  Corrections 

The  following  modifications  to  the  present  beam  segmentation 
method  of  ray  tracing  are  incorporated  in  the  new  ray  tracing  method. 

1 .  New  Ray  Data  Extraction  Process 

In  view  of  the  error  propagating  from  the  Neuman  boundary,  the 
injected  ray  data  plane  (IRD  plane)  was  made  a  user-defined  parameter. 
Previously,  it  was  fixed  to  coincide  with  the  axial  location  of  the 
Neuman  boundary  at  the  end  of  the  simulation.  Now  the  user  specifies 
its  axial  location  depending  on  the  length  of  the  run.  Subsequent 
simulation  runs  8  cm  long  made  to  verify  emittance  growth  showed  the 
solution  to  be  accurate  at  a  distance  equal  to  one  period  of  the  mag¬ 
netic  field  away  from  the  boundary.  This  is  illustrated  by  the  values 
of  Table  1.  The  difference  between  the  tunnel  emittance  values  at  this 
distance  is  less  than  1  percent. 

This  change  in  the  location  of  the  IRD  plane  required  a  corre¬ 
sponding  modification  in  the  logical  diagram  of  Fig.  2.  The  modified 
logical  diagram  is  shown  in  Fig.  26.  Two  ray  data  extraction  G 
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Table  1.  Tunnel  emittance  and  beam  filling  factor  values 
obtained  using  the  new  ray  tracing  method. 


Fig.  26.  Modified  logical  diagram  for  the  new  ray  tracing  method. 


boundaries  are  placed  at  the  axial  location  where  ray  data  will  be 
extracted.  These  G  boundaries  are  configured  according  to  Fig.  27.  The 
potential  values  at  the  nodes  of  these  two  boundaries  will  be  used  to 
start  the  succeeding  continuation  runs.  The  first  G  boundary  marks  the 
position  where  the  next  simulation  run  will  be  started.  The  second  G 
boundary  provides  the  node  potentials  on  the  x-plane  IRD.  These  poten¬ 
tials  are  crucial  in  the  proper  launching  of  the  rays  in  the  continua¬ 
tion  run. 

The  radial  and  axial  coordinates  of  each  node  falling  on  these  G 
boundrles  must  be  exactly  specified  in  the  program,  since  they  will  be 
the  basis  for  the  following  continuation  runs.  All  axial  nodes  are 
separated  by  a  distance  of  Az,  as  shown  in  the  figure.  The  radial 
coordinates  are  obtained  through  the  code.  After  setting  the  number  of 
nodes  that  will  be  used  to  describe  the  height  of  the  tunnel  in  the 
logical  diagram,  the  code  is  used  to  calculate  the  radial  coordinates  by 
specifying  mesh  generation  only  (CONTROL  ■>  ’M’)  and  specifying  a  print¬ 
out  of  the  mesh  generation  result  (PRINT  -  'L')  in  a  trial  simulation 
run. 


2.  New  Ray  Data  Injection  Process 

The  Neuman  boundary  at  the  beginning  of  the  tunnel  in  the  continu¬ 
ation  run  was  replaced  by  a  B  boundary.  The  logical  diagram  is  shown  in 
Fig.  28.  Again,  the  two  G  boundaries  close  to  the  Neuman  boundary  at 
the  end  of  the  tunnel  aire  ray  data  extraction  boundaries.  The  B  bound¬ 


ary  and  the  combination  B/G  boundary  at  the  beginning  of  the  tunnel  are 
ray  data  injection  boundaries.  They  are  configured  according  to  Fig. 


B 


B 
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Fig.  27.  Configuration  of  the  additional  G  boundaries 
in  the  modified  logical  diagram  of  Fig.  26 


I 


Fig.  28.  Modified  logical  diagram  for  the  continuation 
runs  in  the  new  ray  tracing  method 


( 


29.  The  potentials  at  the  ray  data  extraction  boundaries  of  the  pre¬ 


vious  run  are  transposed  here  to  start  the  simulation.  The  hybrid  B/G 
boundary  is  used  in  the  determination  of  space-charge  distribution. 
Since  the  rays  are  not  launched  at  the  beginning  of  the  tunnel,  erron¬ 
eous  potentials  will  be  set  at  the  nodes  falling  on  the  IRD  plane 
because  of  the  absence  of  rays  in  the  region  between  the  B  and  B/G 
boundaries.  Defining  the  nodes  lying  on  the  IRD  plane  to  be  B  nodes 
fixes  the  potential  at  these  nodes  and  bypasses  the  charge  deposition 
process . 

In  order  that  the  same  logical  diagram  could  be  used  for  all 
continuation  runs,  all  simulations  are  started  at  z  -  0.  This  is  not 
only  convenient,  but  it  also  avoids  any  effect  that  a  change  in  the 
logical  diagram  can  bring  in  the  solution.  Again,  the  radial  and  axial 
coordinates  of  the  nodes  lying  on  the  vertical  B,  B/G,  and  two  G  bounda¬ 
ries  must  be  exactly  specified  just  like  in  the  ray  data  extraction 
process . 

3.  Magnetic  Field  Corrections 

The  error  resulting  from  the  magnetic  field  values  is  easily 
corrected  by  performing  linear  interpolation  on  the  measured  B-field. 
This  will  give  us  discrete  zero  magnetic  field  points  that  will  be 
inserted  in  the  B-field  data  set.  Table  2  shows  the  magnetic  field 
values,  with  the  zero  fields  points  incorporated. 
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Table  2.  Magnetic  field  values  of  the  B-field  profile  of  Fig.  8 


z  (m) 


0 .oooooooc^oe 

1 .40«4000C-02 
1.4221000C-02 
1 .434I00OC-02 
1 .447S0O0C-02 
1 .4C02000C-03 
1 . 4724000C-02 
1 .4ISi000e-02 
1 .4*t3000C-02 
l.SllOOOOC-02 
1.52320001-02 
1.51440001-02 
1 .54510001-02 
1 . 54110001-02 
1 . 52450001-02 
1.51220001-02 
1 .55950001-02 
1 .41240001-02 
1.42510001-02 
1 .43400001-02 
1 .45020001-02 
1 .44140001-02 
1.42410001-02 
C-02 
1.20150001-02 
1 .21420001-02 
1 . 22490001-02 
1 .23940001-02 
1.25210001-02 
1.24500001-02 
1.22220001-02 
1 . 29040001-02 
1 .10110001-02 
1-0  2 
1.42450001-02 
1 .44120001-02 
1 .45190001-02 
1-02 
1 .42910001-02 
1 .49200001-02 
1 .90420001-02 
1 .91240001-02 
1.91010001-02 
4240001-02 
1 .95550001-02 
1 .94420001-02 
1-02 
1- 

2.00410001-03 
2 .01900001-03 
2 .01120001-02 


.05210001-02 
.04940001-03 
.04250001-03 
.09520001-03 
. 10290001-02 
. 12040001-02 
. 11110001-02 
.  H40000C-02 


B  (T) 


O.OOOOOOOl^OO 

o.oooooeoi»oo 

2.12119241-0< 
1.32325*51-05 
2.24251*11-05 
1.41242**l-05 
4  .*1431*11-05 
<.05421*01-05 
2.022112*1-0* 
*.111014*1-05 
1.1040*541-04 
1.122*4241-04 
1.551*1*51-04 
1  .*2242141-04 
2.12252511-04 
3.14554221-04 
2.54420111-04 
2.04154201-04 
3.414044*1-04 
3.44452421-04 
4.22*124*1-04 
4. *2*51041-04 
5.55409411-04 
4.512225*1-04 
2.240521*1-04 
*.29*52521-04 
*.*41254*1-04 
1.12*4*2(1-01 
1.11044*41-01 
1.4*022421-01 
1.20(005(:-01 
1.91191291-01 
2.1*01(111-01 
3.19499121-01 
2.20410411-01 
1.0*514221-01 
.422 
.222 
4.(2224221-01 
5.42412241-01 
4.11510*11-01 
2.110120*1-01 
2.9*9124*1-01 
4 .*2041*21-01 
* .(14(1](1-01 
1.011424*1-02 
1.10420421-02 
1 .1242*411-02 
1.2514**11-02 
1.11*1**41-02 
1 .1(451041-02 
2 

1 .51500241-02 
1.54522091-02 
1.42*91(11-02 
1. *4901411-02 
1.21*24141-02 
1.2202*221-02 
1 .42250(11-02 
1 . ((0(2441-02 


Z(m) 

B  (T) 

3.1S*20*0t-*2 

1. *5204221-03 

2.12l4000t-*2 

3.03*20411-*3 

2.1*410001-03 

3.11524211-02 

3.1*4*0001-03 

3.241*25*1-02 

3.30*50001-02 

2.1*4201(1-03 

2.32330001-03 

2.542(3541-03 

3.214*00*1-03 

2.21**4211-03 

2.342(0001-03 

2. *11*4251-02 

2.3*010001-03 

3.1*1**131-02 

2.32100001-03 

3.4**04421-03 

2.24520001-02 

3. *10*5521-02 

2.2*440001-03 

4.20454911-02 

3.11110001-02 

4.(02*0411-02 

3.121(0001-03 

5.34445421-02 

3.11(50001-03 

5.01159221-02 

2.14*20001-03 

*.1144(011-02 

2.141*0001-02 

2.04*21111-02 

2.12440001-03 

2. *1102411-02 

2.1(210001-03 

0.5455**41-02 

2.40000001-02 

*.154*1211-02 

2.41220001-03 

1 .014*92(1-01 

2.42540001-03 

1 .11511151-01 

2.41(10001-03 

1.1*442011-01 

3.450(0001-03 

1.24052(11-01 

2.4(150001-02 

1.2*041*51-01 

3.42(20001-03 

1.12125(51-01 

3.4*4*0001-03 

1.314454(1-01 

3.501(0001-03 

1.12441101-01 

2.51410001-03 

1.39921411-01 

2.52200001-03 

1.24214211-01 

2.51*20001-02 

1.10112111-01 

2.55240001-03 

1.0**52251-01 

3.5(510001-02 

1 .014*9421-01 

2.52240001-02 

9.225(7511-02 

3.59050001-02 

0.504241(1-02 

3.(0120001-02 

2. ((02(421-02 

2.(1590001-03 

*.*4101721-02 

2.(2(40001-03 

4 .04(42(01-02 

3.(4110001-02 

5.2235400C-02 

3.(5400001-02 

4. *4(10151-02 

2.(4(20001-03 

4.114210*1-02 

3.(2*40001-02 

1.(2*01411-02 

2.4*210001-03 

1.11(0(951-02 

3.204*0001-03 

2.4*50(1(1-02 

3.21250001-02 

3.110(1201-02 

2.21020001-03 

1 .*(2*1(01-02 

2.242*0001-03 

1.(555*1(1-02 

2.25540001-03 

1.1241(091-02 

2.24*10001-02 

*.4***9021-01 

2.2(100001-03 

5.02220521-01 

2.29120001-02 

2.912025(1-0] 

2.(0440001-02 

1 .01259921-05 

2.10994001-02 

0 .000000014^00 

2. 11910001-02 

-2.(2*29911-05 

2.111(0001-02 

-4.(45)4241-01 

2 .14450001-02 

-*  .*9009221-01 

2 .15220001-02 

-1.44(12(21-02 

2 .44990001-02 

-2.002(9)21-02 

2.(4240001-02 

-2 .40954951-02 

-1.141712)1-02 

Table  2.  Continued 


Z  (m) 

B  (T) 

Z  (m) 

B  (T) 

2.»OIOOOOC>02 

-3 .•4i4419e-03 

9.(700000C-02 

-1.392a(4aE-01 

2.9207000C-02 

-4 .(001441C>02 

3 .(•27000C-02 

-1.4a97235E-01 

2.9334000C-02 

-5.4191423C-02 

3.(994000C-02 

-1  .(2(9593e-01 

2.»4(1000C-02 

-(.293S99SC-02 

9.7091000E-02 

-1.74727900-01 

2 .95IIOOOC-a2 

-7.3290219C>02 

3.7309000C-02 

-1  .aaaaoo(o-oi 

2.971S000C-e2 

-• .41444(tC-02 

9.7339000B-02 

-2  .0ia(972C-01 

2.9I42000C-02 

-9.S907<<4C-02 

3.74(20009-02 

-2 .12241(1E-01 

2 .99<9000C-02 

-1 .0I91494C>01 

3.75t9000C-02 

-2.2105741E-01 

3 .009(000C-02 

-1 .237S3(9C-01 

3.771(OOOC-02 

-2.2a9a341E-01 

3 .0223000C-02 

-1.394C034C>01 

3.7a43000B-02 

-2 .33(5(99E-01 

3 .03SOOOOC-02 

-l.S9011t0C-01 

3 .7970000C-02 

-2.3(724210-01 

3.0473000C-02 

-1.703St90C-01 

3.I097000C-02 

-2.3751200E-01 

3 .OC04000C-02 

-l.t409904C-01 

3  .I224000C-02 

-2.3(04SaiE-01 

3.0731000C-03 

-1.97(S«01C-01 

9.I391000C-02 

-2.32a(920E-01 

3 .OISIOOOC-02 

-2.0t(7407e-01 

3.I47IOOOC-02 

-2.2(732230-01 

3 .09ISOO0C-O2 

-2.1t31244C-01 

3.I(05000B-02 

-2.ia31344B-01 

3.1113000e-02 

-2.2(73223e-01 

3.I732000C-02 

-2.0a(7407B-01 

3 .1239000e-03 

-2.32l<920e-01 

3.II99000C-02 

-1.97(5(010-01 

3 . 13C<O0Oe-O2 

-2.3<04S91C-01 

3.191(0001-02 

-1.94099040-01 

3.1493000e-02 

-2.3791200C-01 

3 .9113000C-02 

-1.70399900-01 

3 . lC20000e-02 

-2.3C72421C-01 

3 .9240000C-02 

-1 .5501iaOB-01 

3 . 1747000e-02 

-2.334SC99C-01 

3.93(7000B-02 

-1.394(0340-01 

3 .lt74000e-02 

-2.3t9t341t-01 

3 .9494000B-02 

-1 .21791(9E-01 

3 . 2001000C-02 

-3.210S741C-01 

3 .9(21000C-02 

-1.09914940-01 

3 .212l000e-02 

-2.12241ilC-01 

3 .974IO0OB-O2 

-9.9507(940-02 

3.22SS000C<02 

-2 .01IC972C>01 

3.9I75000B-02 

-9.41444(90-02 

3.23I2000C-02 

-l.llll00<c>01 

4.0002000C-02 

-7.32502190-02 

3 . 2J09000C>02 

•1.74727I0C>01 

4.0129000C-02 

-(.29399990-02 

3 . 3C3(000C-02 

«1.<2(9993C-01 

4 .029(000t-02 

-5.41914220-02 

3.27C3000C-02 

-1.4I97239C-01 

4 .03I3000C-02 

-4.(0014410-02 

3.2l90000e-02 

-1.3S29C4IC-01 

4.0910000C-02 

-3.94(44190-02 

3 . 30l7000e-02 

•1.2190229C-01 

4.0(37000C-02 

-3.1(37173E-02 

3 .3144000C-02 

-1.0912149C-01 

4.07(40000-02 

-2.(0954990-02 

3 . 3271000C>02 

•9.S(3401(E-02 

4 .OI91OO0C-O2 

-2.00799370-02 

3.339I00OC-O2 

•  t  .449214IC-02 

4.101I000E>02 

-1.44(17970-02 

3.392$OOOC>02 

-7.4<12303C-02 

4.1145000E-02 

-».(900927E-03 

3 . 3C J2000C-02 

-«.«l75354C-02 

4 .1272000E-02 

-4.(453424E-03 

3 . 3779000C-02 

-S. 91210921-02 

4.1399000E-02 

-2.(2979910-05 

3 . 390(000C-02 

-S.0999199C-02 

4.14(2500E-02 

0.00000000»00 

3.4033000C-02 

-4 .49l9777t-02 

4.192(000E-O2 

2.(2979910-05 

3 . 41(0OOOC-O2 

-4 .019Sll<C-02 

4  .1(53000E-02 

4.(4534240-03 

3 . 42I7000C-02 

-3 .7102S32C-02 

9.17I0OOOE-O2 

9.(9009270-03 

3.4414000C-02 

-3 .440S992C-02 

4 .1907000E-02 

1.44(17970-02 

3.4941000C-02 

-3.2097<37C-02 

4 .2034000E-02 

2.00799370-02 

3.466I000C-O2 

-3 .0717942C-02 

4 .21(1000E-02 

2.(0954950-02 

3 .4799000C-02 

-2 .990131IC-02 

4.22II0OOE-O2 

3.1(371730-02 

3.4922000C-02 

-3 .0717942C-02 

4.2419000E-02 

3.94(44190-02 

3. J049000C-03 

-3 .2097(17C-02 

4.2942000E-02 

4.(0014410-02 

3 . 91 7i000C-02 

-3 .4409I92C-02 

4.2((9000E-02 

5.41914220-02 

3 . 5303000C>02 

-3 .7102932C-02 

4.279(000E-02 

(  .2915999E-02 

3 .5430000C-02 

-4.01991t6C-02 

4 . 2923000E-02 

7.1250219E-02 

3 . 5597000C-02 

-4 .49I9777C-02 

4 . 3090000E-02 

9.41444(90-02 

3 . 56I4000C-02 

-S.0999199C-02 

4.3177000E-02 

9 .5507((4E-02 

3.5I11000C-02 

-5.I121092C-02 

4.3304000B-02 

1.0a91494E-01 

-<.(179394C-02 

4.3431000B-02 

1.23753(90-01 

4.399I000B-02 

1.394(034E-01 

■•!■"!  51 

4  .  1(I9000E-02 

1 .5501ia0E-01 

-9 .9<340l4E-02 

4 . 3ai2000B-02 

1.7035990E-01 

-1.0912149C-01 

4.  1939000E-02 

KWTFTTrnsn 

3.(573000C>92 

-1 .2190229C-01 

4 . 40((OOOE-02 
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Table  2.  Continued 


Z(m) 


4 .411 3000C>02 
4.4)20000C-02 
4.4447000C>02 
4.4S74000C-02 
4 .4701000C>02 
4  .4I2I000C>02 
4.4fSS000C-02 
4  .S042000C-02 
4  .S209000C-02 
4 .Sl)<000C-02 
4  .S4(3000C-02 
4.SSt0000C-02 
4  .S717000C-02 
4  .SI44000C-02 
4.Sf7l000C*02 
4 .(094000C-02 
4 .(32S000C-02 
4  .<3S2000C-02 
4  .<479000C-02 
4 .(<0<000C-02 
4  .<733000C-02 
4 .<I(0000C*02 
4  .<9I7000C>02 
4 .7114000C-02 
4 .7241000C*02 
4.73<IO0OC-O3 
4.7495000C*02 
4 .7<23000C-02 
4 .7749000C-02 
4 .7l7<000C-02 
4  .•003000C-03 
4.I130000C-02 
4  .•257000C-02 
4  .I3I4000C-02 
4 .•511000C-02 
4  .•C3I000C>02 
4  .•745000C-02 
4 . II92000C-02 
4 .9019000C-02 
4 .914C000C-02 
4 .9373000C>02 
4 .9400000C-02 
4 .9937000C-02 
4 .9<S4000C-03 
4 .9741000C-02 
4 .9904000C>02 
S.003S000C-02 
».01C3000C>02 
3.03I9000C-03 
3.041(000C-03 
3.0343000C-02 
3.0<70000C-03 
3.0797000C-03 
3.0934000C-e2 
3. 1031000C>e3 
3.117I000C>92 
3.1303000C>02 
3.1433000C-03 
3. 1359000C>03 
3.1(l<000C-03 


B  (T) 


2.0I(7407C>01 
2.1I31344C-01 
2.2(73323C-01 
3.33IC930C>01 
2.3(04341C-0l 
2.3731200C-01 
3.3<72421t-01 
2.33(3(99C-01 
3.2I3I341C-01 
2.2103741C-01 
2.12241<1C-01 
2.01I<972C>01 
1 .IIII00<C-01 
1 .74727t0C-01 
1 .<2(9333C-01 
1 .4I97233C-01 
1 .332I<4IC-01 
1 .2130223C-01 
1 .0912143C-01 
9.3<3401(C>02 
2 

7.4<12103C>02 

<.<173334C-02 

S.I121032C-02 

S.0333199C-02 

4 .49I9777C-02 

4 .01931I(C*02 

3.7102332C-02 

3.44034922-02 

3.20974372-02 

3.07179422-02 

2.99013112-02 

3.07179422-02 

3.20974372-02 

3.44031922-02 

3.71023322-02 

4.01931142-02 

4 .49497772-02 

3.03331992-02 

5.41210322-02 

4.41753342-02 

7.44121032-02 

4.44921442-02 

9.54340142-02 

1.09121432-01 

3.21302232-01 

3.33244442-01 

3.44972332-01 

3.42493332-01 

1.74727402-01 

3.44440042-01 

2.01449722-01 

2.12241412-01 

2.21037412-01 

2.24513412-01 

2.33434992-01 

2.34724212-01 

2.37312002-01 

2.34043112-01 

2.32149202-01 


Z  (m) 


11130002-02 

39400002-02 

20470002-02 

21940002-02 

23210002-02 

24440002-02 

23730002-02 

27020002-02 

24290002-02 

29340002-02 

30130002-02 

32100002-02 

33370002-02 

34440002-02 

33910002-02 

37110002-02 

34430002-02 

39720002-02 

40990002-02 

42240002-02 

43330002-02 

44000002-02 

44070002-02 

47340002-02 

47975002-02 

41410002-02 

49440002-02 


B(T) 


24732232 

30332442 

00474072 

97454012 

44099042 

70354902 

55011402 

39440342 

23753492 

00914542 

55074442 

43444412 

32502132 

29359952 

41914222 

40034412 


34373732-0 

40954932-0 

00709372-0 

44417172-0 

49009372-0 

44334242-0 

42479912-0 

00000002«0 

42179912-0 

44534242-0 


4.  R95  Value  Corrections 


The  error  caused  by  the  erroneous  R95  values  is  avoided  by  recal¬ 
culating  the  R95  values  in  the  continuation  runs.  The  correct  values 
are  obtained  by  using  the  beam  current  present  in  the  initial  simulation 
in  all  computations. 

B.  Simulation  Using  the  New  Method 

These  modifications  and  corrections  brought  about  a  stair-step 
method  of  beam  segmentation.  In  this  new  approach,  the  segments 
included  in  each  simulation  run  are  made  to  overlap  each  other.  Figure 
30  illustrates  this  overlapping  process.  Simulation  'a*  covers  the 
segment  from  the  cathode  up  to  a  distance  of  approximately  8  cm  down¬ 
stream.  Simulation  'b*  (starting  and  ending  at  z  *  5.5  cm  and  z  -=  13.6 
cm,  respectively)  overlaps  and  resolves  the  end  section  of  simulation 
•a’,  allowing  us  to  discard  the  solution  adjacent  to  the  Neuman  bound¬ 
ary. 

Simulation  performed  in  Section  III  was  repeated  using  this  new 
beam  segmentation  method.  The  resulting  beam  profile  is  shown  in  Fig. 
31a.  It  is  remarkably  in  good  agreement  with  the  beam  profile  of  Fig. 
12.  The  plots  of  the  beam  filling  factor  and  tunnel  emittance  for  the 
beam  profile  of  Figs.  12  and  35  (shown  in  Figs.  32  and  33)  clearly 
illustrate  this  agreement  between  the  two  solutions.  To  further  illus¬ 
trate  the  improved  accuracy  resulting  from  this  new  method  of  beam 
segmentation,  the  curves  of  Figs.  32  and  33  are  plotted  along  with  the 
beam  filling  factor  and  tunnel  emittance  obtained  using  the  old  beam 
segmentation  method.  These  curves  are  shown  in  Figs.  3ii  and  35. 
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Fig.  31.  Electron  trajectories  using  the  new  method  to  verify  tunnel  emlttance  growth 


Fig.  32.  Beam  filling  factor  versus  axial  distance  using  (a)  the  new  method 
with  beam  segmentation,  and  (b)  simulation  without  segmentation. 


Fig.  33.  Tunnel  emlttance  versus  axial  distance  using  (a)  the  new  method 
with  beam  segmentation,  and  (b)  simulation  without  segjmentation 


Tunnel  emittance  versus  axial  distance  using  (a)  the  new  method  with  beam  segmentation 
(b)  simulation  without  segmentation,  and  (c)  the  old  method  with  beam  segmentation. 


Fig.  35.  Beam  filling  factor  versus  axial  distance  using  (a)  the  new  method  with  beam  segmentation 
(b)  simulation  without  segmentation,  and  (c)  the  old  method  with  beam  segmentation. 


VI.  THE  NEW  RAY  TRACING  METHOD 


This  section  is  a  formal  presentation  of  the  ray  tracing  method 
that  incorporates  the  modifications  discussed  in  the  previous  section. 
The  step-by-step  procedure  is  demonstrated  through  an  example.  Brief 
discussions  are  included  to  clarify  and  bring  up  important  details  that 
will  allow  DEMEOS  users  to  perform  ray  tracing  without  any  difficulty. 

A.  Ray  Tracing  Procedure 

The  previous  section  presented  a  simulation  using  this  new  method 
of  beam  segmentation.  That  simulation  will  be  referred  to  as  each  step 
of  the  ray  tracing  procedure  is  presented.  Figures  enclosed  in  paren¬ 
theses  are  actual  values  used  in  the  simulations. 

The  ray  tracing  procedure  is  divided  into  two  steps: 

1.  Modeling  and  simulation  of  the  gun  and  its  magnetics. 

2.  Extraction  and  injection  of  ray  data  in  the  continuation  runs. 
Since  the  difference  between  the  new  and  the  old  method  of  ray 

tracing  lies  mainly  in  the  continuation  process,  only  the  process  of 
extracting  and  injecting  of  ray  data  will  be  presented.  Modifications 
that  will  affect  the  simulations  in  the  gun  region  will  be  pointed  out 
as  they  come  up. 


1 .  Ray  Data  Extraction 

1.  Adjust  the  simulation  length  to  include  at  least  four  periods  of  the 
magnetic  field.  In  the  example,  the  first  single-peak  and  four 
periods  of  the  double-peak  magnetic  fields  were  included  in  the 
first  simulation  run.  This  is  shown  as  segment  'a'  in  Fig.  30. 
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2.  Define  x-plane  IRD  to  be  the  axial  coordinate  of  the  next  to  last 
zero  B-field  location.  (This  is  at  z  -  5.^17975  cm  in  the  example.) 
Note  that  at  least  one  extra  field  point  after  the  last  zero  B  field 
is  needed  to  properly  run  DEMEOS.  Also,  all  zero  B-field  locations 
must  be  discretely  defined  in  the  magnetic  field  data  set. 

3.  Modify  the  logical  diagram  of  the  gun  region  by  adding  two  G  bounda¬ 
ries  in  the  tunnel  region,  as  shown  in  Fig.  26.  The  G  boundary 
closest  to  the  Neuman  boundary  at  the  end  of  the  region  should 
coincide  with  the  location  of  the  x-plane  IRD  (the  plane  where 
injected  ray  data  will  be  extracted).  The  other  G  boundary  is 
positioned  at  a  distance  of  Az  away  from  the  x-plane  IRD,  as  shown 
in  Fig.  27.  This  distance  is  calculated  as  follows: 

Az  -  L/N 


where 


L  -  length  of  simulation  run 

N  »  node  number  defining  the  axial  end  of  the  run 

Note  that  the  recessed  nodes  (odd  numbered  J’s)  are  positioned  at  a 
distance  of  Az/2  from  x-plane  IRD. 

In  the  example,  node  No.  215  was  set  for  the  end  of  the  first 
region  at  z  «  8.1^675  cm.  Dividing  this  value  of  z  by  215  gives 
a  Az  of  0.03789  cm.  Dividing  the  axial  distance  of  the  x-plane  IRD 
(z  »  5.^7975  cm)  by  Az  gives  li|^.62.  Rounding  off  gives  node  No. 
1^45  for  both  the  x-plane  IRD  and  the  second  G  boundary.  With  this 
value  of  Az ,  the  second  G  boundary  has  axial  coordinates  of  5.^7975 
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era  for  even  numbered  J’s  and  5.^9870  cm  for  odd  numbered  ones. 
Similarly,  the  first  G  boundary  (at  node  No.  has  axial  coor¬ 

dinates  5. ^^186  cm  and  5.^6080  cm,  respectively. 

Tighten  the  tolerance  to  increase  the  accuracy  of  the  solution.  The 
following  convergence  tolerances  were  used  in  the  example: 

ALPHA-1  -  1.500  ALPHA-2  -  L.50  ALPHA-2  -  1.300 

EPS-1  -  0.003  EPS-2  -  0.01  EPS-V  =  0.003 

Note  that  convergence  to  within  these  specified  tolerances  may  not 
be  completely  met.  In  simulation  runs  made  for  this  study,  it  was 
observed  that  data  away  from  the  end  of  the  run  will  reach  stability 
even  if  the  tolerances  are  not  completely  satisfied.  Simulation 
could  be  terminated  once  the  solution  past  the  x-plane  IRD  becomes 
stable  and  solutions  up  to  the  x-plane  IRD  can  be  considered  to  have 
converged.  It  will  be  helpful  to  define  at  least  two  x-planes  after 
the  x-plane  IRD  so  that  comparison  of  the  solutions  in  the  overlap¬ 
ping  segments  could  be  made.  Figure  30  illustrates  how  the  simula¬ 
tion  runs  overlap.  The  end  segments  of  each  run  are  included  in 
succeeding  simulation  runs. 

2.  Ray  Data  Injection 

1.  Adjust  the  simulation  length  to  include  at  least  four  periods  of  the 
magnetic  field.  In  the  example,  six  periods  of  the  double-peak 
periods  are  included  in  the  first  continuation  run,  resulting  in  a 
simulation  length  of  8  cm.  This  is  shown  in  Fig.  30  as  segment  'b*. 
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2.  Outline  the  tunnel  according  to  the  logical  diagram  of  Fig.  28.  The 
two  G  boundaries  near  the  Neuman  boundary  follow  the  discussion  in 
the  data  extraction  process. 

3.  Start  the  simulation  at  z  -  0.  This  means  that  the  vertical  B 
boundary  will  have  an  axial  coordinate  of  a  -  0  cm  for  even  numbered 
J's,  and  Az/2  (0.03789  cm/2)  for  the  odd  numbered  J's. 

4.  Configure  the  B/G  boundary  according  to  Fig.  29.  The  B  and  G  nodes 
have  axial  coordinates  equal  to  Az  (0.03789  cm)  and  1.5  Az  (0.05684 
cm),  respectively. 

5.  Subtract  the  beam  voltage  potential  (26,000  V)  from  the  potentials 
at  each  of  the  nodes  on  the  two  G  boundaries  discussed  in  the  data 
extraction  process.  Transpose  the  resulting  values  to  the  nodes  of 
the  B  and  B/G  boundary.  The  first  G  boundary  corresponds  to  the 
vertical  B  boundary,  while  the  second  G  boundary  corresponds  to  the 
combination  B/G  boundary. 

6.  Since  the  rays  are  not  launched  at  the  start  of  the  simulation,  the 
restriction  that  the  magnetic  field  be  zero  at  the  start  of  the 
simulation  does  not  apply.  However,  the  field  must  be  zero  at  the 
location  where  the  rays  are  launched  (at  z  -  0.03789  cm).  At  least 
two  extra  field  points  (one  before  the  start  of  the  region,  and 
another  after  the  end)  are  needed  to  properly  run  DEMEOS. 

7.  Tighten  the  tolerance  to  Increase  the  accuracy  of  the  solution.  Set 
the  convergence  tolerances  according  to  the  values  given  in  step  4 
of  the  ray  data  extraction  process. 

8.  Reset  the  x  values  of  the  rays  in  the  injected  ray  data  file  to 
zero. 


-  66  - 


Since  we  are  simulating  segments  of  considerable  length,  it  will 
take  several  iterations  before  the  solution  reaches  convergence.  Again, 
we  do  not  have  to  attain  complete  convergence  of  the  entire  solution. 
As  soon  as  the  solution  past  the  x-plane  IRD  becomes  stable,  convergence 
of  the  solution  can  be  assumed.  This  is  one  aspect  of  this  beam  segmen¬ 
tation  process  that  is  very  advantageous.  We  are  not  interested  in  the 
solution  past  the  IRD  plane,  since  they  are  going  to  be  discarded  any¬ 
way.  Again,  it  is  recommended  that  at  least  two  x-planes  after  the  x- 
plane  IRD  be  defined  to  help  us  in  determining  the  stability  of  the 
solution. 

To  reach  a  convergence  in  the  potential  of  0.003  percent 
(Av.  ./v.  .  <  0.00003)  sometimes  requires  an  abnormal  amount  of  computer 

^  iJ  1  iJ 

time.  To  avoid  excessive  run  times,  it  is  recommended  that  LOOPMAX  be 
set  equal  to  a  value  less  than  or  equal  to  10  and  set  SAV  ■  'B*.  This 
saves  the  output  of  the  simulation  after  the  10th  iteration,  allowing 
evaluation  of  intermediate  solution.  To  start  the  11th  loop/iteration, 
set  SAVD  ■  'Y'.  This  Instructs  the  code  to  use  the  result  of  the  last 
iteration  as  initial  condition.  Remember  that  an  output  and  input  file 
must  be  defined  in  these  processes. 

All  continuation  runs  are  set  up  in  similar  fashion.  No  changes 
in  the  logical  diagram  is  needed  as  long  as  the  segments  included  in 
each  simulation  are  kept  the  same. 
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VII.  CONCLUSIONS 


As  mentioned  in  the  beginning  of  this  report,  the  need  to  verify 
emittance  growth  is  the  motivation  behind  this  study.  Because  of  the 
inconsistencies  discovered  in  the  present  ray  tracing  method,  highly 
accurate  simulation  over  the  full  length  of  the  tube  was  not  possible. 
This  thwarted  our  investigation  of  beam  behavior  and  compelled  us  to 
deviate  from  our  initial  objective.  We  proceeded  to  investigate  the 
causes  of  the  inconsistencies  in  the  ray  tracing  approach  and  developed 
a  method  that  could  be  used  in  verifying  a  very  sensitive  parameter, 
such  as  tunnel  emittance.  Now  that  we  are  able  to  accurately  trace  the 
beam  over  the  full  length  of  the  tube,  we  could  proceed  to  verify  occur¬ 
rence  of  emittance  growth. 

A .  Verification  of  Emittance  Growth 

Simulation  of  the  entire  length  of  the  tube  was  made  by  continuing 
the  simulations  performed  in  Section  V.B.  The  beam  was  traced  up  to  a 
distance  of  approximately  30  cm  downstream.  The  beam  profile  resulting 
from  this  simulation  is  shown  in  Figs.  36a,  b,  c,  d,  and  e.  The  length 
of  the  segments  and  the  periods  of  magnetic  field  included  in  each 
simulation  run  are  depicted  in  Fig.  30.  Figures  37  and  38  show  the 
plots  of  the  beam  size  and  tunnel  emittance  for  this  beam.  Tables  3  and 
show  the  values  used  in  these  two  figures.  It  also  demonstrates  that 
the  distance  of  1  period  of  magnetic  field  away  from  the  Neuman  boundary 
is  sufficient  to  obtain  error  free  solutions. 
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Fig.  36.  Electron  trajectories  using  the  new  method  to  verify  tunnel  emlttance  growth. 


(c)  Second  continuation 
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(e)  Fourth  continuation 
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Fig.  37.  Beam  filling  factor  versus  axial  distance  for  beam 


Table  3.  R95/A  values  obtained  from  the  beam  of  Fig.  36. 


Z(cm) 

R95/A 

(a) 

(b) 

(c) 

(d) 

(e) 

2.1 

0.53012 

3.1 

0.6l76i< 

i<.l 

0.38i<85 

5.1 

0.68209 

6.1 

0.H8H30 

0.i<8iJ9'< 

7.1 

0.H9350 

0.i<9939 

8.1 

0.62570 

0.5761J5 

9.1 

0.58652 

10.1 

0.i<6522 

11.1 

0.73263 

0.732^8 

12.1 

0.J<801U 

0.480'<i4 

13.1 

0.56792 

0.56086 

U.l 

0. 632214 

15.1 

0.57719 

16.1 

0.449879 

17.1 

0.60281 

0.600744 

18.1 

0.448126 

0.448393 

19.1 

0.59511 

20.1 

0.526449 

21.1 

0.446100 

22.1 

0.52809 

0.52801 

23.1 

0.595443 

0.59591 

2H.1 

0.4414422 

0.440302 

25.1 

0.58199 

26.1 

0.449863 

27.1 

0.55821 

28.1 

0.61933 

29.1 

0.52938 

Table  4.  T-einittance  values  obtained  from  the  beam  of  Fig.  36 


Z(cm) 

1 

■-Emlttance 

(a) 

(b) 

(c) 

(d) 

(e) 

2.1 

0.01633 

3.1 

0.01240 

4.1 

0.01445 

5.1 

0.01750 

6.1 

0.01620 

0.01597 

7.1 

0.01620 

0.01621 

8.1 

0.02060 

0.01848 

9.1 

0.01714 

10.1 

0.01876 

11.1 

0.02319 

0.02320 

12.1 

0.01932 

0.01928 

13.1 

0.01992 

0.01969 

14.1 

0.02256 

15.1 

0.02024 

16.1 

0.01963 

17.1 

0.01881 

0.01896 

18.1 

0.01566 

0.01753 

19.1 

0.01880 

20.1 

0.01817 

21.1 

0.01494 

22.1 

0.01815 

0.01811 

23.1 

0.01523 

0.01514 

24.1 

0.01558 

0.01530 

25.1 

0.01439 

26.1 

0.01836 

27.1 

0.01961 

28.1 

0.02013 

29.1 

0.01788 

The  plot  of  the  beam  filling  factor  over  the  full  length  of  the 
tube  shows  the  average  size  of  the  beam  remaining  fairly  constant. 
Tunnel  emittance  is  shown  increasing  from  a  minimum  value  of  0.012  to  a 
maximum  value  of  0.023.  Since  this  was  a  nonthermal  simulation,  correc¬ 
tions  will  still  have  to  be  made  to  be  able  to  compare  these  values  to 
that  predicted  by  Eq.  8  in  reference  2. 

Figure  37  also  shows  tunnel  emittance  exhibiting  a  wavelike 
behavior.  After  growing  to  a  maximum  of  about  0.023,  tunnel  emittance 
proceeded  to  decay  to  a  value  of  approximately  0.01^.  It  then  started 
growing  again  until  it  reached  a  value  of  approximately  0.021  at  the  end 
of  the  tube.  If  tunnel  emittance  is  indeed  growing  and  decaying  in  a 
cycle,  this  could  be  a  very  significant  discovery.  This  implies  that 
tubes  could  be  designed  such  that  the  dc  beam  will  have  minimum  tunnel 
emittance  at  the  end  of  the  tube.  This  will  result  in  an  additional 
Increase  in  beam  transmission  and  high  tube  efficiency. 

It  is  noteworthy  to  point  out  that  tunnel  emittance  did  not  decay 
to  its  minimum  value.  After  reaching  its  first  maximum,  it  decayed  to  a 
value  of  approximately  0.01^.  Even  if  it  is  truly  cycling,  it  seems 
that  it  does  not  return  back  to  its  initial  value. 

B.  Suggestions  for  Further  Study 

We  have  now  predicted  emittance  growth  over  the  full  length  of  a 
tube.  The  simulation  that  was  made  resulted  in  several  interesting 
results  that  need  experimental  verification.  The  next  step  before 
embarking  on  any  experimental  measurements  is  to  extend  the  simulation 
until  determination  of  the  wave-like  behavior  could  be  made.  Also, 
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perturbations  in  the  entrance  condition  could  be  introduced  to  determine 
how  these  will  affect  tunnel  emittance.  In  particular,  the  position  of 
the  beam  waist  with  respect  to  the  entrance  magnetic  field  could  be 
varied  to  determine  if  there  is  an  optimum  position  that  will  result  in 
minimum  emittance  growth  and  higher  beam  transmission.  The  effects  of 
shunting  the  magnetic  field  could  also  be  simulated  by  varying  magnetic 
field  values  along  the  tube.  Adjustments  of  the  B-field  could  be  made 
at  the  locations  of  the  beam  excursions  in  Fig.  36  to  observe  the 
effects  of  this  process. 

To  determine  the  accuracy  of  this  new  beam  segmentation  method, 
the  5637  gun  and  the  5638  body  could  be  put  together  and  tested  in  the 
gun  analyzer.  A  thermal  simulation  using  this  new  method  has  to  be  made 
so  that  accurate  comparison  between  measured  and  predicted  values  could 
be  made.  Once  accuracy  and  reliability  of  this  new  method  is  estab¬ 
lished,  a  tube  could  be  built  to  verify  experimentally  the  predicted 
emittance  growth. 

C.  Final  Remarks 

In  addition  to  the  development  of  a  new  method  of  ray  tracing  for 
emittance  growth  studies  in  very  long  beams,  the  work  uncovered  what 
appears  to  be  a  fundamental  problem  in  all  codes  which  use  a  paraxial 
approximation  for  space-charge  (see  Section  IV. B).  The  implication  is 
that  codes  employing  such  paraxial  approximation  will  have  inherent 
errors  in  the  solution.  This  can  be  verified  by  comparing  the  results 
of  paraxial  based  codes  to  nonparaxial  results  from  DEMEOS. 
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